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Motivated by a possibility to optimize modelling of the population evo¬ 
lution we postulate a generalization of the well-know logistic map. Gener¬ 
alized difference equation reads: 

Xn+i=rxl{l-xl), ( 1 ) 

X € [0,1], (p, q) > 0 , n = 0 , 1, 2,..., where the two new parameters p and 
q may assume any positive values. The standard logistic map thus corre¬ 
sponds to the case p = q = 1 . For such a generalized equation we illustrate 
the character of the transition from regularity to chaos as a function of r 
for the whole spectrum of p and q parameters. As an example we consider 
the case for p = 1 and q = 2 both in the periodic and chaotic regime. We 
focus on the character of the corresponding bifurcation sequence and on 
the quantitative nature of the resulting attractor as well as its universal 
attribute (Feigenbaum constant). 


1. Introduction 

Since the dawn time the chaos is an indispensable part of human life. 
One of the most famous models with chaotic behavior is the logistic map. 
Its simple analytic form meant that it was used in many scientific disciplines 
such as biology, cryptography, communication chemical physics and stock 
market [DEliail!. One of the mathematical models of chaos was discovered 
by Mitchell Feigenbaum (1978) O [6l [71 [8] . Feigenbaum considered ordinary 
difference equations used for example in biology to describe the development 
of population in its dependence on time. He discovered that population 
oscillates in time between stable powers (fixed points), the number of which 
doubles according to changes in the power of external parameter. The 
generalization of the Feigenbaum model contains all first-order difference 
equations f{xn) = Xn+i- The condition for the existence of chaos is the 
single maximum of a function f{xn)- Feigenbaum also proved that the 
transition to chaos is described by two universal constans a and 5, later 


( 1 ) 
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named the Feigenbaum numbers. 

We quote here a few definitions which are basic in the development of this 
exposition. 

Definition 1. [cf. [9], p. 62] 

Let I = [a, b] C a < b, f{I) C I. Discrete iterations of a function 
f : I ^ I are called functions : I ^ I defined inductively: = id, 

r+' = /or,n = 1,2,3,.... 

Definition 2. [[S], p. 62] 

The sequence {xn}, where Xn = f"'{xo) is called the sequence of iterates of 
a function / generated by point xq G I. 

Definition 3. [cf. [9], p. 65] 

A point s £ I, with /(s) = s is called a first-order fixed point of /. A fixed 
point s of a function / is locally stable if there exists a neighbourhood Us 
such that all sequences of iterates are convergent to s for xq G Us. 
Theoreml. [see [9], p. 68] 

If function / is differentiable in fixed point s and ]/ (s)| < 1, then s is 
attractive. 

Definition 4. [cf. [9], p. 72] 

Let f : I ^ I. We say that s is a periodic point of the function /, or that 
s generates a cycle of order fc > 2, if f^{s) = s and / s. 

Theorem2. [TO] 

If function / has a cycle of order different to 2^, then / is chaotic. 
Theorems, [cf. [9], p. 79] 

For arbitrary continuous function /:/—)•/ existence cycles of order m 
involve the existence cycles of order n in the following order: 

3 » 5 » 7 » 2 • 3 » 2 • 5 » 2 • 7... » 2^ • 3 » 2^ • 7 » ... » 
2^ >> ... >> 8 >> 4 >> 2 >> 1. 

The conventional logistic map 

The Logistic map 

f{Xn) = Xn +1 = rXn{l “ Xn), X £ [0,1], {p, q) > 0, n = 0,1, 2,... (2) 

is one of the most simple forms of a chaotic process. Basically, this map, like 
any one-dimensional map, is a rule for getting a number from a number. 
The parameter r is hxed, but if one studies the map for different values 
of r (up to 4, else the unit interval is no longer invariant) it is found that 
r is the catalyst for chaos. After many iterations x reaches some values 
independent of its starting value 3 regimes: r < I : x = 0 for large n, 
I < r < 3 : X = constant for large n, 3 < r < Too : cyclic behavior, where 
Too ~ 3.56, Too < r < 4 : mostly chaotic. 

An interesting thing happens if a value of r greater than 3 is chosen. The 
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map becomes unstable and we get a pitchfork bifurcation with two stable 
orbits of period two corresponding to the two stable fixed points of the sec¬ 
ond iteration of /. With r slightly bigger than 3.54, the population will 
oscillate between 8 values, then 16, 32, etc. The lengths of the parame¬ 
ter intervals which yield the same number of oscillations decrease rapidly; 
the ratio between the lengths of two successive such bifurcation intervals 
approaches the Feigenbaum constant 5 = 4.669.... The period doubling 
bifurcations come faster and faster ( 8, 16, 32, ...), then suddenly break off. 
Beyond a certain point, known as the accumulation point r^o, periodicity 
gives way to chaos. In the middle of the complexity, a window suddenly 
appears with a regular period like 3 or 7 as a result of mode locking. The 
3-period bifurcation occurs at r = 1 -|- 2\/2, and period doubling then begin 
again with cycles of 6,12,... and 7,14,28,... and then once again break off to 
chaos. 


2. Generalized logistic map 

One-dimensional representations are the simplest arrangements capa¬ 
ble to producing a chaotic movement. Let us consider p,q > 0. We postulate 
non-linear function: 


fr{x) = rx^{l — x^), X G [0,1] 
set by difference equation: 


( 3 ) 


Xn+i = rxl{l - xl), X G [0,1], n = 0,1,2,.... (4) 


For the first time the function ([3]) has been presented in HU. The max- 

/ \i/'? / 

imum of the function ([3]) is reached for x = ( j , because (x) = 
r{px'P~^ — {p + q)x'P~^’^~^) and //(x) = 0 x = ^ > 3^ 9 (P; O') > 0. 

It results that and for r G fo, ^ ^ ^^ j 





(p+g) 


= Traax function fr is /([O, 1]) C [0, 1]. 


Parameter Vmax depending on p and q accepts the following limit values 
(FigHD: 

lim rmax = 1, lim r^ax = oo 

q^ao q^0+ 

lim rmax = OO, lim r^ax = 1 

p-^oo p-^0+ 
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Fig. 1: Values of r^ax for generalized logistic map depending on p and q 
(left) and depending on k (right). 


For the sake of parameters p and q the function (3) generates a lot of 
new functions which can be analyzed using mathematical methods, jnst like 
for a logistic map, thus obtaining very interesting results. 

a) Let us assume that k = - and k > 0. Then r is depend on the value k 

fc+i 

and r G 0, 

b) If A: = 1, then p = q and rmax is independent on p, q and rmax = 4 
always - for logistic map, we have p = q = 1 and rmax = 4. 

k + l 

If A: —>■ oo, then rmax 1, because lim^^oo = 1- 


c) 


d) 

e) 


If A; —)• 0+, then r„ 


oo, becanse lim^^oo 


k + l 

(fc + 1)^ 


= + 00 . 


Now let us consider a case, when p = 1 and q > 0. Then we have 

fr{x) = rx{l — x'^), xG[ 0, 1], n = 0,1,2,.... (5) 


and r G 


2±i 

p, (g+1) 

n 


Taking into account the above relationship we can observe change max¬ 
imum of the function (5) depending on q (Figl2]). 
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A) The property of Xn+i = rxn(l — Xn) for p = 1, q = 2 in the 
periodic region. 

In the following part of this paper we will show that the Feigenbaum 
model of transition to chaos is correct for representation (4) for given pos¬ 
itive values of parameters p and q. Therefore, we will study the dynamics 
of this map in the case of p = 1 and q = 2. Then the representation (5) 
assumes the following form: 

fr{x) = rx{l — x'^), a:€[0,1], n = 0,1,2,... (6) 

A function fr in the point x = ^ has a local maximum and the value of 
function in this point is equal to /r(^) = ■ From this relationship, we 

have /([0,1]) C [0,1] for r € (0,^^]. The trajectories behavior generated 
by fr depends on r and, therefore, we have to distinguish several intervals 
of this value: 

1) For r G (0,1] all the sequences of iterates xi,X2,X3, ...,Xn depend on 
the initial value xq tend to s = 0. This point is the fixed point of the map 
for r G (0,1] and according to Definition 3 it is also the stable point of this 
map. 

2) Let r = ri > 1. If x = rx{l — x'^) then, xi = 0, X 2 = \J so for r > 1 

function (6) has two fixed points. Using Theorem 1 we get: (x) = r—3rx^, 

I//(0)1 = r and \ fr{^J^)\ = |3 - 2r|. 

Therefore, we have that xi = 0 is an unstable fixed point, whereas X 2 = 
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is an attractive fixed point for r E ( 1 , 2 ). 

Change r from 1 to 2 has caused the migration of the attractive point from 


0 to 


C2 


Conclusion 1. 

For 1 < r < 2 there exists a stable fixed point attracting all points. When 
r = ri = 2, then \ fr = |3 — 2r| = 1 and this point stops being an 

attractive point. 

3) For r = r 2 > 2 we will look into the stability of fixed points fr and fr 
of the map (5) as a function of parameter r, because the Feigenbaum model 
is generated by bifurcations connected with iterations of a function. 

Let us consider the second iteration of fr, namely fr‘^\x) = r'^x — r'^x^ — 

+ 3r^x® — Sr^x"^ + with fixed points: xi = 0 , X 2 = ^3 = 


4 

r 


xa = \l\ + Because fr — 3r^x^ — Sr'^x'^ + 

15r'^x‘^ — 21r^x® + 9r'^x® we have fr^‘^\xi) = r^, fr^‘^\x 2 ) = (3 — 2r)^, 
fr^'^\x 3 ) = fr^'^\xA) = 9 — 2r^ whereas for 2 < r 2 < V5 fixed points X 3 
and X 4 are attractive points according to condition in Theorem 1 . 
Conclusion 2. 

The fixed point x = of is also the fixed point of as well as of 

all higher iterations. 

Conclusion 3. 

If the fixed point of fr^'^ becomes unstable, then it is also the unstable fixed 
( 2 ) 

point of fr ^ and of all next iterations. 

From the inequality |/'(s)| > 1 we have \f'^'^\s)\ = \f' [f{s)]f' {s)\ = 

I/'(^)I>1- 

If \/5 < rs < r 4 , then the fixed points of fP'^ become repulsive at the same 
time. Following this instability, the fourth iteration shows two new pitchfork 
bifurcations giving cycle 2 ^ order for four attractive fixed points, which is 
called period duplication. 

To generalize the above examples we get: 
a) For r„_i < r < there exists a stable cycle 2”“^, whose elements 
xo, xi, X 2 , X 3 ,..., X 2 n-i _4 are defined by the following connections: 


frixi) = Xi+I, fP \xi) = Xi, lYifriXi)] < 1 

i 


b) For Tn all points of the cycle of 2” ^ become unstable simultane¬ 
ously and the bifurcation creates a new stable cycle of 2 ” order for fP = 
fr~\fr~P existing for r„_i < r < Xn- 
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c) For a map (6) the value of parameter Vn+i we can describe the following 
connection: r„+i Ri -^3 + where tq = 1, ri = 2 ,r 2 = \/5, •■•,5’oo- 

d) There exists a point of accumulation of an infinite number of the bi¬ 

furcation of period duplication for the finite value r which we appointed as 
roo- roo = lim^^oo'Tn = ~ 2.303. 

Taking into account the above, Vn+i ~ Vn, \/3 + r = r, — r — 3 = 0, 
thus r = 

e) For the representation (6) there exists a constant, whose value corre¬ 
sponds to the Feigenbaum constant: 

6 = lim R.4.6I. 

n^QO rn+1 - Tn 

Proof. Let us assume that r^-i = a. Then = -^3 -|- Vn-i = \/3 + a, 

rn+i = V3 + Tn = \/3 + v^3"+a and 6 = lim^^oo = \/l4 -k 2 a/13 « 

4.61. 

A table of 2” cycles type and values of is given below. It is clearly vis¬ 
ible that doubling bifurcations come faster and faster (8, 16, 32, ...) and 
suddenly break off beyond a certain point known as the accumulation point 
and periodicity gives way to chaos. 


n 

1 

2 

3 

4 

5 


11 

12 


00 

2” cycle 

2 

4 

8 

16 

32 


2048 

4096 


accum. point 

Pn 

2 

2.23606 

2.28824 

2.29961 

2.30209 


2.3027755 

2.30277562 


2.30277563... 


Table 1: The algebraic orders of the values of for n = 1,2,... are given 
by 1,2,4,8,.... 


B) The property of Xn+i = rxn(l — Xn) for p = 1, q = 2 in chaotic 
region. 

a) For Too the bifurcation sequence ends with a set of the infinitely nu¬ 
merous points, which is called the set of Feigenbaum attraction. 

b) If Too < r < then we observe the irregular dynamics of map (6) 
with narrow ranges in which the set of attraction has a periodical character. 

c) Numeric calculations for the representation (6) show that the largest 
range is for r = Tc ~ 2.451. Then there exists 3-cycle (Figi3]), which accord¬ 
ing to th.2 and th.3, implicates the chaos as well as the existence of cycles 
of different orders. 

d) For the map (6) at r = rmax = 


3v^ 


^n(l ^n) — frmaxi^n), 


^n+1 — 


2 


( 7 ) 
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Fig. 3: Iterates of map (6) strting from xq = 0.99. Left: In the stable 3-cycle 
region for r = Xc, Right: In the chaotic region for r = Vc — Ar, Ar = 0.001 


can be solved by the simple change of variables: 


Xn = cos (yn) = giVn)- 


Then (7) can be converted into 


cos(y„+i) = ^^cos(yn)[l-cos^(yn)] = cos (y„) sin^(y„), 


which has one solution yn+i = arccos cos (y„) sin^(i/ji) and 

-yo + yZ-l + 2/0^ 


( 8 ) 


Xn = arccos (- 


'/ 3 {-yo + \/-l + yo^ 


T + 


V3 


The invariant density function p^{x) of fs^jx) can be calculated using 

“ 2 “ ~T~ 

the following dehnition: 


^ N-l ^ N -1 

lim — > 6 (x — Xn) = lim — > 6 (x — qiyn)) 
N^OO N ^ ^ ’ N^oo N ^ ^ y\yn)) 

n=0 n=0 


Using p{y) = 1, (9) becomes: p^{x) = dyp{y)5[x - h{y)] 



—- n 

100 


(9) 
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i.e. 


/ ,_ \ 2/3 

f \/x2 — 1 — x\ — 1 

= I 

^ —T — x>—(y/x'^ — 1— x) — ^ —^273 +1 

y \ / 1—xj ' 

( 10 ) 

The above function is shown in FigSJ 


p(x) 




0.1 0.3 0.5 0.7 0.91 


Fig. 4: The invariant density function of map (6) for r = rmax = Left: 
An analytical form represented by the equation (10). Right: Numerical 
simulation - the histogram for 50000 iterations. 


The invariant density function of the logistic map for r = 
given by . | . 

ttW {l—x)x 


= 4 is 


C) Ordered chaos - construction of an attractor for 
Xn+i = rxn(l - x^) and r = r^ax = 

At the end of this considerations we will present a method of the order¬ 
ing of chaos. A lot of numeric experiments have led to the situation where 
it is possible to rank chaos in such a way that the successive iterations gen¬ 
erated by chaotic dynamic system would reach an attractor of a regular 
shape. This method also has it, that from the sequence of several thousand 
iterations (for representation = rxn{4 — x^) where r = rmax = 

defined recursively) we draw up the following points {xn+i,XnXn+i) on a 
plane - we construct an attractor. In this way an attractor comes into be¬ 
ing, which is visited irregularly by points of successive iterations (FigJS]). 
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Fig. 5: The set of attraction for representation (6) [60000 iterations]: a) for 
f = ?’max(i'ed), r = 2.56(blue), r = 2.53(green), r = 2.45(orange); b) for 
r = rmax(red), r = 2.42(black), r = 2.4(blue). 


D) Bifurcation diagram and Lyapunov exponent. 


Nowadays the Lyapunov exponent is applicable in various scientihc dis¬ 
ciplines such as biology, engineering, bioengineering and informatics [laiig 
[iziiisiiigEn]. A quantitative measure of the sensitive dependence to initial 
conditions is the Lyapunov exponent, which is a measure of the exponential 
separation of nearby orbits. For discrete time system Xn+i = f{xn), for an 
orbit starting with xq the Lyapunov exponent A(xo) can be calculated by 


A(xo) = lim - ^ln|/'(xi)|, (11) 

i=0 

A positive Lyapunov exponent can be considered as an indicator of chaos, 
whereas negative exponents are associated with regular and periodic be¬ 
havior of the system. Figure [6] shows the bifurcation diagram and Lya¬ 
punov exponent for the equation ([6|). It is clear that even in chaotic areas 
(^'cxD < r < there are many periodic intervals (A < 0). 
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r 


Fig. 6: The bifurcation diagram and the Lyapunov exponent for the map 
© and 2.2 < r < 


3. Summary 

Currently, modeling real-world systems is very popular. A well-known 
model that takes into account the deterministic chaos is the logistic map. 
However, this one-dimensional map can be controlled only by one parameter. 
In the present contribution, we postulate a generalization of the classical 
logistic map which is generally controlled by three parameters. This could 
allow to adapt to the conditions prevailing in the modeled system. 

We focus on the Vmax parameter (depending on p and q) which is re¬ 
sponsible for the dynamics of the system - among others shown that '\ip = q 
then r,nax = 4. 

As an example, we carry out analytical and quantitative analysis in the 
case where p = 1 and q = 2 both in the periodic and chaotic regime where 
the Lyapunov exponent has positive values. It turns out that the dynamics 
of this equation is faster than in the case of the logistic map - the value of 
the parameter r^ax ~ 2.3 is less than the value of r^ax = 4 for the logistic 
map. In the periodic regime we identified Feigenbaum constant 6 which is 
typical for all dissipative systems (also for the logistic map). For chaotic 
area we have found an analytical form of the invariant density function. For 
both regimes we have proposed specific form of data representation which 
allowed to obtain a non-trivial structure of an attractor set. 

Footnotes 

This work was partially supported by the Centre for Innovation and Transfer 
of Natural Sciences and Engineering Knowledge (University of Rzeszow). 
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